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We report an intrinsic form of the inverse spin Hall effect (ISHE) in ferromagnetic (FM) metal with Rashba 
spin orbit coupling (RSOC), which is driven by a normal charge current. Unlike the conventional form, the 
ISHE can be induced without the need for spin current injection from an external source. Our theoretical 
results show that Hall voltage is generated when the FM moment is perpendicular to the ferromagnetic layer. 
The polarity of the Hall voltage is reversed upon switching the FM moment to the opposite direction, thus 
promising a useful readback mechanism for memory or logic applications. 



I. INTRODUCTION 

Recent research showed that Rashba spin orbit cou- 
pling (RSOC) at the surfaces of metals can be en- 
hanced by the presence of heavy atoms^'^ and/or sur- 
face oxidation^ in adjacent layers. This inter facial en- 
hancement enables significant RSOC effect to be man- 
ifested in ferromagnetic metals with small or moderate 
atomic number at room temperature^, whereas previ- 
ously, strong RSOC effect is confined only to semiconduc- 
tor heterostructures. In this paper, we choose a typical 
ferromagnetic (FM) metal (Co) as the central conducting 
layer, sandwiched between an oxide and a Pt layer, the 
latter supplying the heavy atoms [see Fig. 1(a)]. The 
electron accumulation which develops in the FM layer in 
the presence of a charge (unpolarized) current is theoret- 
ically evaluated via the non-equilibrium Green's function 
(NEGF) method in the ballistic limit. In the presence of 
s-d coupling, the incoming charge current becomes polar- 
ized by the FM moments in the central region. When this 
intrinsic polarization of current is coupled to the RSOC, 
an inverse spin Hall effect (ISHE)^'^ will be induced. 
Thus, a Hall voltage is generated without the need for 
spin injection from an external spin polarizing layer. By 
contrast, in previous works, the ISHE is experimentally 
realized by injecting spin polarized current^'^ from an 
external FM electrode, or by the inflow of pure spin 
current^'^^^, generated externally e.g. via spin pumping 
or non-local spin accumulation. In this work, we show 
theoretically that Hall voltage can be generated when 
the FM moment in the central region is oriented perpen- 
dicular to the plane, which persists at room temperature. 
Furthermore, the generated Hall voltage can be reversed 
symmetrically when the FM moment is switched to the 
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Electrode Edge width 

FIG. 1. (a) Schematic diagram of the proposed FM mo- 
ment detector utilizing the ISHE phenomenon, (b) Lattice 
discretization of the device for tight-binding NEGF calcula- 
tion (top view) . The Hall voltage Vt is the potential difference 
between the two strip electrodes running along the top and 
bottom edges of the central region. The covered area schemat- 
ically shows the edge width which will be used for the Hall 
voltage calculation. 



opposite direction. Thus, the charge current-induced 
ISHE signal can be used to detect the polarity of the FM 
moment, and potentially serve as a read-back mechanism 
in memory applications. 



II. MODEL HAMILTONIAN AND THEORY 

The schematic diagram of the FM moment detec- 
tor is shown in Fig. 1(a); the central region com- 
prises of a triple-layer structure for the enhancement of 
RSOC within the FM (Co) layer^. Charge accumula- 
tion within the Co layer is calculated via tight-binding 
NEGF method^ ^. To perform the tight-binding calcula- 
tion, the central region of the device is discretized into 
(M X TV) lattice of points [see Fig. 1(b)]. The conduction 
electrons within the Co layer experiences the RSOC ef- 
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feet and the s-d exehange interaetion with the local FM 
moments M(^,0). Thus, the Hamiltonian of the central 
region can be expressed as He = Hk-\-Hm-\-Hrso^ where 
Hk is the kinetic term, Hm the s-d coupling term, and 
Hrso the RSOC term. The Hamiltonian can be expressed 
as": 

Hx = ^ ^ [^tdl^^^dmna 
mna 

— t{dl^_^i^^^dmna + <^ln,n+lcr^^ncr + h.C.)], (1) 

Hm =Y1 sgn[cr]Mcos(i9)d]^^^d 

mna 

mna 

+Msm{e)e'-Ul^,d^na, (2) 

Hrso = —itso[{dl^_^i^^dmn — <^ln-l,n^wn) ^ O'y 

mnaa' 

~i^m,n-\-l^'mn — <^L,n-l^^n) ^ CTx], (3) 
I 



where M denotes the s-d coupHng strength, tso = 
denotes the RSOC strength. Similarly, the Hamiltonian 
of the normal metal (NM) leads, and the coupling energy 
between the leads and central region can be expressed as: 

^L{R) = ^ [^'^(AnnaO^rnna 
mna 

~ ^i^m-\-l,na^'mna + ^m,n-\-la^rnna + h.C.)]. (4) 

= "^[tWona^lna + tR^M+l^na^^Mna + h.C.].(5) 
na 

From the eigenvalue equation of the total Hamiltonian 
and the definition of retarded Green's function, one can 
obtain an equation: {E — Hmn + ^^)^m^(^^) = ^• 
From this relation, one can obtain a series of linear equa- 
tions involving G]^^^{aa) by considering each spatial 
point im^n). For instance, within the central region, i.e. 
1 < m < M, one obtains: 



I=[E-At- aMcos{e)]G^^^^{aa) - e^^'^ M sm{e)G^^^^{aa) + t[G2-i,^{<y<y) + G::;+i,„.(aa) 



Collectively, all these equations can be expressed in ma- 
trix form: {E[I] - [H])[GY = L The infinitely large 
matrix [H] consists of sub-matrices denoting the Hamil- 
tonian of the central region {[He]) and the couphng co- 
efficients between the two leads and the central region 
separately {[tl/r^c])- Following standard procedures in 
the tight-binding method, one then obtains: 

iE[I]-[Hc]-[ni-m[Gcr = I, (7) 
in which the non-zero terms of the self energy are: 

\yn,n -\r _ M^^^ In^]^'^' M^''^' 

[^m,m\L{R) ~ r Jm,0(M+l) J 0(M+1) ,0(M+1) r Jo(M+l) ,m * 

The retarded Green's functions of the isolated left (right) 
lead [^^]o(7Vf^i) o(m+i) expressed as (for the left 

lead): 
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■^^Xi(Pn)e*"Xi(Pn')- 



(8) 



In the above, ki is the wave vector along the semi-infinite 
longitudinal direction, Xi{Pn) is the ith eigenfunction in 
the transverse dimension at site (0,n) in the lead, which 
can be expressed as: 



Xi{Pn) 



1 2 


iTTTl 




- sin . 



N+l N+1 



(9) 



The retarded Green's function of the central region can 
then be solved by inverting Eq. (7). Thus one can 
express the lesser Green's function [G]^ via the Lan- 
greth formula: [G]< = [G]^[E]<[G]^ in which [E]< = 
E/.=L,i?(P]M ~ with being the Fermi distri- 

bution function within lead /i. The total charge accumu- 
lation for a given cell at lattice coordinate (m, n) with an 



area of is given by: 



"2^ 



Tr[G]<(E)dE, (10) 



where the trace is over the spin degree of freedom. The 
surface charge density pm,n is then given by pm,n = 

The Hall voltage at longitudinal position x = m (Vtm) 
is given, up to a proportionality constant, by the differ- 
ence in the surface charge density between the top and 
bottom edges corresponding to x = m, i.e. Vtm ^ ^Pm = 
Ptm — Pbm- In calculating the charge densities ptm and 
pbrm we consider a finite width of each edge. The val- 
ues of ptm and pbm are averaged over some number of 
rows (W) adjacent to the top and bottom edges, such 
that Wa ^ 0.3 nm. Thus, the surface charge density 
difference along x-direction is given by: 



Apr, 



,N-n-\-l 



w 



(11) 



In practice, the Hall voltage Vt is given by the potential 
difference between the two electrodes. We assume that 
each electrode runs along the entire length of the edges 
(i.e., from m = 1 to M). Thus, computationally, the Hall 
voltage Vt is given by the surface charge density difference 
averaged over the longitudinal dimension, i.e.. 



Apa 



M 



1 2^n=l 



~ Pn 



M xW 



(12) 
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FIG. 2. (a) Distribution of AyOm along the longitudinal x- 
axis. The FM moments in the central region are oriented 
along the zLx (blue), ±2/ (red) and ±2; (black) directions. The 
corresponding electron density distributions are schematically 
depicted in (b), (c) and (d), respectively. Fig (e) shows a pos- 
sible electrode configuration to detect the Hall voltage when 
FM moment is along ±x directions. The following parameter 
values are assumed: RSOC strength of tso = 3 eV, bias volt- 
age of y = 4 eV, and s-d coupling strength of M = 0.85 eV. 
The central region is discretized into a lattice of (200 x 100) 
unit cells. 



III. RESULTS AND DISCUSSION 

The following parameters are assuraed in the numerical 
calculations: (i) The device is modeled at room temper- 
ature (T = 300 K); the Fermi energy of the central FM 
layer is set to 7.38 eV, which is a typical value for Co^^. 
(ii) The lattice cell dimension is set to a = 0.045 nm, 
which is significantly smaller than the Fermi wavelength 
(a ^ A/10), so that the lattice Green's function model 
can simulate a continuum system to a good approxima- 
tion, (iii) The coupling strength is t = ^^^2 = 18.69 
eV, while for simplicity, the coupling between the lead 
and the central region is set to t^/ji = O.St, (iv) The 
RSOC strength in Eq. (3) is given by tgo = f^- For 
a typical FM RSOC material, a lies between 4 x 10~^^ 
and 3 x 10~^^ eVm^'^^'^^, which translates to a range 
of coupling parameter values of 0.4 < tgo < 3.3 eV. (v) 
The s-d exchange energy is set to \M\ = 0.85 eV^^. (vi) 
The electrochemical potentials of the two leads are set 
to fiL = —fJ^R = 2 eV. (vii) Finally, the central region is 
discretized into a square lattice of (M x TV) = (200 x 100) 
of unit cells. This corresponds to an actual dimensions 
of (9 nmx 4.5 nm) for the central region. 

We first calculate the transverse charge density dif- 
ference of Apm as a function of the longitudinal posi- 
tion X, when the FM moments in the Co layer are sep- 
arately oriented along ±x, ±y and ±z directions. The 
results are plotted in Figure 2(a). Our results show 
that when the FM moments are in the ^/-direction, the 
spatial charge distribution pm,n is symmetric about the 
central longitudinal axis of n = (TV + l)/2, resulting 
in Apm = 0. When the FM moments are switched to 
—y direction, the charge distribution pm,n remains sym- 
metric about the central longitudinal axis, and hence 
Apm is still [see Fig. 2(b) for a schematic represen- 
tation]. Thus, when the FM moments of the Co layer are 
aligned along ±y, no Hall voltage would be observed. By 
contrast, when the FM moments are oriented along the 
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FIG. 3. The spatial distribution of pm,n (in unit of ^). 
The longitudinal and transverse dimensions are expressed in 
unit of the lattice constant a. The s-d coupling strength is 
M = 0.85 eV, the RSOC strength is tso = 3 eV, the bias 
voltage is y = 4 eV, the central region is a (200 x 100) lattice 
(9 nm x4.5 nm). 
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FIG. 4. The detailed distribution of ptm and pbm (in unit of 
^) at (a) the top, and (b) bottom edges of Fig. 3. The 
longitudinal and transverse dimensions are expressed in unit 
of the lattice constant a. 

x-direction, Apm is symmetric about the central point 
(m,n) = ((M + l)/2, (A/" + l)/2) [see Figs. 2(a) and 
(c)]. Furthermore, the sign of Apm is reversed when the 
FM moments are switched to the —x direction. However, 
when averaged along the entire edge, the surface charge 
density difference will be Apav = due to the point sym- 
metry. However, if the Hall electrodes extend to only half 
the entire length of the top and bottom edges [as shown 
schematically in Fig. 2(e)], a finite Hall voltage can still 
be detected. 

Of greater interest is the case where the FM moments 
are along the out-of-plane z-direction. The charge den- 
sity difference Apm is symmetric about the central ver- 
tical axis, i.e. m = (M + l)/2. Thus, a finite Apav, i-e., 
a Hall voltage Vt is generated [see Figs. 2(a) and (d)]. 
Since only charge current but not spin current is injected 
into the system, the above can be regarded as a charge 
current-induced ISHE in FM metal with RSOC. 

In the following, we will investigate this charge current- 
driven ISHE in greater detail. The spatial distribution 
of 

Pm,n is plotted over the central region when the FM 
moments are in the +2;-direction [see Fig. 3]. For clarity, 
the detailed distribution of charge densities ptm and pi^m 
are shown in Figs. 4(a) and (b). It is observed that the 
surface charge density is larger along the top edge, which 
will result in a finite Hall voltage or ISHE effect. 

Fig. 5(a) shows the dependence of Apav on the RSOC 
strength tso for different edge width W . Apav, and hence 
the Hall voltage Vt across the central region, show an in- 
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FIG. 5. (a) The oscillatory increase of /S.pav with the RSOC 
strength tso, for different edge width W and FM moment 
orientations (solid lines for dashed lines for —z ). The 
central region is discretized into a (200 x 100) lattice, with 
dimensions (9 nmx4.5 nm). The s-d coupling strength is 
M = 0.85 eV, while the bias voltage is y = 4 V. (b) The 
schematic diagram of the spin electron distribution due to 
Yang- Mills- like Lorentz force. The number of Hall deflection 
pairs increases with increasing RSOC strength. 



to the —z direction. The resulting large difference in the 
Hall voltage corresponding to the two FM orientations 
(±z) suggests that the ISHE can be utilized in a metal 
FM RSOC system for the sensitive detection of the FM 
moment orientation. 

In summary, we have investigated the inverse spin Hall 
effect (ISHE) which is induced by the combination of 
RSOC and s-d coupling to the FM moments. A Hall 
voltage is generated when the FM moments are oriented 
in the perpendicular-to-plane direction. The Hall volt- 
age increases in an oscillating manner with the RSOC 
strength tso- The polarity of the Hall voltage is reversed 
when the FM moment is switched to the opposite direc- 
tion. This property suggests the utility of the ISHE in 
FM metals with strong RSOC effect for the detection of 
the FM moment direction, e.g., as a possible memory 
readback mechanism. 



creasing trend with the RSOC strength tso-, but in an os- 
cillatory manner. The physics underlying this oscillatory 
increase can be understood in terms of the Yang-Mills- 
like Lorentz force arising from the RSOC gauge^^'^^. 
Here, the FM moments M merely play the role of sus- 
taining a vertical spin polarization of current but do not 
contribute directly to the Lorentz force. The Lorentz 
force leads to the transverse separation of electrons of 
opposite spins [shown schematically in Fig. 5(b)], which 
we term as "Hall deflection pair". Since there are un- 
equal number of charges on the two transverse sides, each 
Hall deflection pair will contribute to a charge Hall volt- 
age. For a fixed M, an increase in the RSOC strength 
results in an increasing number of Hall deflection pairs 
along the length of the device, as shown in Fig. 5(b). 
The charge imbalance in a Hall deflection pair coupled 
with the increase in the number of such pairs with the 
RSOC strength provide a heuristic explanation of the os- 
cillatory increase of Ap^^ and the Hall voltage with the 
RSOC strength. The sign of Ap^v can be reversed by 
switching the orientation of FM moment between ±z. 
Furthermore, Ap^v is not sensitive to the definition of 
edge, i.e. the general trend remains unchanged for the 
range of edge width W considered in our calculation. 

Previous work on RSOC in semiconductors has shown 
that when the surface charge density difference is in the 
order of lO^^e/m^, the generated Hall voltage will be 
large enough for detection (0.1 mV)^^. The electron den- 
sity in our metallic FM RSOC device will be much higher 
than that in semiconductors, so that Ap^^ could attain 
a value of the order of lO^^e/m^ and generate a sizable 
Hall voltage oiVt ~ IV . By selecting optimal tso which 
corresponds to the peak Hall voltage values (see Fig. 5), 
we conjecture that a reasonably large Apav^ hence Vt can 
be measured when the FM moments are oriented along 
the out-of-plane z axis. The charge density difference 
Apav switches in sign upon reversal of the FM moments 
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